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Abstract. We present a systematic search for Fe Ka emission from young stellar 
objects of the p Ophiuchi star forming region observed in the Deep Rho Ophiuchi XMM- 
Newton Observation. 



1. The Deep Rho Ophiuchi XMM-Newton Observation (DROXO) 

The Deep Rho Ophiuchi XMM-Newton Observation, DROXO, is an XMM-Newton 
Large Project with a nominal exposure time of 500 ksec centered on core F of the p Oph 
star forming region. A catalog of 1 1 1 X-ray source s, more than 50 % of which iden- 
tified with known cloud members, was presented by iPillitteri et all (120101) . That paper 
also comprises a detailed description of the general data analysis steps including data 
filtering, source detection, treatment of the background etc. Here we present the re- 
sults of a systematic search for Fe Ka emission in all 111 DROXO X-ray sources. The 
Fe Ka emission of one DROXO source, the Class I p rotostar Elias 29, has already been 
discussed in a separate paper dGiardino et all 12007 ). Our analysis also includes this 



object and we confirm the results by Giardino et al. 



2. Fe Ka emission in Young Stellar Objects 



The X-ray spectra of Young Stellar Objects (YSOs) can be explained by thermal emis- 
sion from their coronae that are a few times 10 6 K hot. In roughly a dozen YSOs an 
emission line at 6.4 ke V has been identified that can not be attribu ted to highly ionized 
coronal material (e.g. Imanishi et al.ll200H : iTsujimoto et alj|2005l) . In the Sun this line 



represents Ka fluorescence of neutral iron in the photosphere irradiated with X-rays 
from the solar corona during flares. In YSOs this line is usually attributed to fluo- 
rescence from neutral iron in the accretion disk following X-ray illumination from the 
star's corona. In some cases the line equi valent widths have been shown to be incompat - 
ible with fluorescence calculations (e.g. lGiardino et aUl2007l : ICzesla & Schmidl2007l) . 
An alternative excitation mechanism for the line is electron impact ionization. Giardino 
et al. suggest the accretion streams between the disk and the stellar surface as possible 
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Figure 1 . Observed Fe Ka equivalent widt hs and X-ray temperature compared to 
calculations for photospheric fluorescence by lDrake et al.l (l2008h . Solid lines repre- 
sent models for different heights, R, of the exciting source in units of stellar radii. 
Error bars for the data points are 68 % confidence levels. 
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location for this process. Another possible scenario that may explain the large observed 
equivalent widths in these cases without the need to invoke collisional excitation is the 
(parti al) occupation of the illuminating continuum flux, e.g. a stellar flare behind the 
limb (DrakeetaL||2008). The identification of the origin of the Fe Ka line in YSOs can 
greatly benefit from time-resolved spectroscopy that enables a direct evaluation of the 
continuum visible at the moment when the line is produced. 



3. Fe Ka spectroscopy of DROXO X-ray sources 



We have carried out systematic spectral fitting of all 111 DROXO sources with XSPEC 
v.12. Spectra were generated for a variety of good- time-intervals, spectral binnings, 
all EPIC detectors individually as well as joint pn+MOS 1+MOS 2, such that for each 
X-ray source a total of 24 time-averaged X-ray spectra are available. Time-resolved 
spectroscopy is performed in intervals defined using a maximum likeliho od technique 



that d ivides the lightcurve in time-intervals of constant intensity (see e.g. IStelzer et al. 



2007). 



We limit the energy range for spectral fitting to 5 - 9 keV. At these high energies 
the photo-absorption, generally an impo rtant factor for YSO s in p Oph (median Nu for 
DROXO sources is 2.3 • 10 22 cm -2 ; see iPillitteri et al.ll2010h . has negligible influence. 
We adapt a one-temperature APEC model plus Gaussian feature with fixed energy of 
6.4 keV to the data. The 6.4 keV feature is considered detected if the normalization of 
the Gaussian is non-zero with 90 % confidence. 
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Figure 2. Equivalent widths of 6.4 keV emission detected in time-resolved spectra 
vs. the time-averaged value. Evidently, the time-averaged equivalent width is a 
single value for a given star. In the figure, for better visibility of all time-resolved 
measurements for a given star, it has been attributed small offsets from its observed 
value. 
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3.1. Results from time-averaged spectra 

The time-averaged spectra include data from 6 satellite orbits. The 6.4 keV line is 
detected in 9 of 1 1 1 DROXO sources: 4 Class I, 4 Class II, 1 Class III objects. 

Searching for the origin of the line, we compar e the observed l ine equivalent 



widths to calculations for photospheric fluorescence (IDrake et all 120081) . The fluo- 
rescence model consists in a single coronal flare illuminating the photosphere from 
a height, R. The flux in the fluorescence line depends essentially on four parameters: 
the iron abundance in the fluorescence medium, the temper ature of the io nizing coronal 
plasma (kT), its height (R), and the heliocentric angle (fl); iBaskol dl979l) . The models 



shown in Fig. Q] are based on the assumption of 9 = 0, i.e. a flare erupting on the star in 
the direction of the line-of-sight, the case that gives the highest fluorescence flux. The 
calculations for R = apply also to a fluorescent disk illuminated by a stellar flare. 

The sample of equivalent widths observed in DROXO is complemented in Fig. Q] 
by YSOs with Fe Ka detections fro m the Chandra Orion Ultradeep Project (COUP) 
discussed by lTsujimoto et al] d2005l ). Most observed equivalent widths are compatible 
with the predictions of the photospheric fluorescence model with a flare at moderate 
height above the photosphere or a flare irradiating a disk. We detect the 6.4 keV line 
from one Class III object, GY 463, where material related to a disk is not a viable option 
for the line emission site (see also Sect i3.2l) . 

3.2. Results from time-resolved spectra 

There are several reasons to expect the Fe Ka flux to be variable in time. Generally, any 
variation of the ionizing spectrum (in terms of temperature or emission measure) has 
an impact on the strength of the line. Specifically, if flares are responsible for its ex- 
citation a correlation between the presence of Fe Ka emission and flaring is predicted. 
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Figure 3. X-ray spectrum with FeKa line fixed at an energy of 6.4 keV for the 
Class III source GY463 representing the time-interval with a large flare. Red is 
the coronal one-temperature model, and green the Gaussian representing the Fe Ka 
emission. Line equivalent width and 68 % uncertainties are written on the top of the 
figure at the line position which is marked with a vertical line. 
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In practice, the relation betw een flares and Ka e mission may be concealed by occul- 
tation effects as proposed by iDrake et ail (120081) . and time-delays between these two 
phenomena may yield geometrical constraints. 

Our time-resolved spectroscopy of DROXO sources results in Fe Ka detections 
in individual time segments for 5 YSOs. In Fig. |2] we show the equivalent width of 
the 6.4 keV feature in individual time segments vs. the equivalent width in the time- 
averaged spectrum for those objects. This comparison demonstrates how mixing dif- 
ferent activity states by averaging over long time intervals yields a smaller equivalent 
width than the actual maximum value achieved during the observation. 

In Fig. |2] we distinguish three activity states (base level, elevated level and flare) 
defined by different intensity in the lightcurve. Remarkably, the 6.4 keV line is detected 
not only during flares but in some ca ses during quie s cence ('base level'). Especially, 
for Elias29 we recover the result of iGiardino et al.l (120071) . Elias29 is also the star 
with the highest equivalent width (up to 700 A). It is the only object in our sample 
that is incompatible with the fluorescence model taking into account the uncertainties 
(Source 38 in Fig.Q]). The second star in which the FeKo- emission is observed outside 
evident flares is IRS 43. This source is discussed in more detail in Sect. 13.31 For the 
only Class III object with Fe Ka emission, GY463, the time-resolved spectral analysis 
shows that the detection of this line is related to a large flare identified in the lightcurve. 
The spectrum corresponding to the time-interval where the line was detected is shown 
in Fig. [3] 
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Figure 4. Three different spectra 
around the 6.4 keV region for IRS 43: 
(top right) time-averaged with fixed 
line position, (bottom left) time- 
averaged with free line position, and 
(bottom right) restricted to the time- 
interval during which the line was 
detected in the time-resolved spec- 
troscopy; see Sect. I3.3I for details and 
caption of Fig. [3] for an explanation of 
the features shown in the figure. 
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3.3. The puzzling Fe Ka emission of IRS 43 



The Fe Ka emission is located at 6.4 keV for neutral and low ionized iron with ver y little 
dependence of the fluorescence energy on the charge state dKallman et al.ll2004l) . This 
justifies our choice of fixing the line position in the spectral fitting process. However, 
in the time-averaged spectrum of IRS 43 (alias YLW 15) this approach results in a clear 
residual between the model and the data that stems from an excess flux at energies 
between the Gaussian (fixed at 6.4 keV) and the coronal iron line at 6.7 keV; see Fig. [4] 
top right. 

IRS 43 is a binary with 0.5" separation. The northern component is identified only 
at radio and mid-IR wavele ngths, i.e. it is probably either a ve ry young protostar or a 
knot in the outflow (e.g. Girar t et all 120001 : lHaisch et all 12002 ). The X-ray emission 
is most likely produced on the brighter object, IRS 43 S, a Clas s I pro tostar according 
to the (unresolved) spectral energy distribution. iGreene & Ladal (120021) detected stellar 
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absorption features in near-IR spectra of IRS 43 that point to a spectral type of K5 IV/V. 
From the near-IR photometry they derive an extinction of Ay ~ 40 ma g and a stellar 
lumin osity of ~ 3L Q leaving 70% of its bolometric luminosity (10 L G ; Wilkin g~et all 
1989) to the accretion process. However, the conclusion on the accretion luminosity 
relies on the interpretation of the second radio and mid-IR source as a jet and not as an 
embedded binary companion that would contribute to the luminosity derived from the 
spectral energy distribution. 

A fit of the DROXO spectrum with free line position results in a line energy of 
6.48 ±0.05 keV where the uncertainty is for a 68 % confidence level (Fig.|4]bottom left). 
This line is apparently produced by significantly ionized ions (FeXIX - Fe XXIII). 
Those ions are dominating the ionization balance of a collisional plasma in a small 
temperature range of 8 - 14 MK. It is therefore difficult to attribute this line to fluo- 
rescence from a photosphere or an accretion disk where such high temperatures are not 
reached. It seems like in this case the fluorescence takes place in the corona itself. How- 
ever, the coronal column density is probably too small to produce the observed equiv- 
alent width of ~ 150 eV. Another hypothesis is that this fluorescence line is emitted in 
denser plasma such as the one present in the chromosphere or transition region. The 
emitting material may then be photo-ionized rather than collisionally ionized, and this 
would allow for lower temperature. Comparison of the observe d line energy with the 



high-resolution profiles for a photo-ionized plasma calculated by lKallman et al.1 (120041) 



suggests an ionization parameter of log £ ~ 2. For this value of £ , the temperature 
of the photoionization calculations is log T [K] ~ 5.2 and the Fe Ka line emissivity is 
dominated by Fe XIX. 

Inspection of the lightcurve and the time-resolved spectroscopy for IRS 43 reveals 
that the Fe Ka line is detected during 4 time-segments. Only one of them shows the 
energy shift. This time-interval corresponds to a small flare during the last orbit of the 
observation. The spectrum of this time segment is shown in the bottom right panel 
of Fig. |4] One of the remaining three time-intervals with Fe Ka emission has a larger 
flare, while the other two represent the quiescent X-ray state. A possible scenario for 
the shifted line energy may be fluorescence from the chromosphere or transition region 
illuminated by a flare occuring near or partially behind the limb. 
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